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1. Introduction 

The purpose of this note is to apply the recent results on semi-classical trace formula , and 
on quantum BirkhofF normal forms for semi-classical Fourier Operators to inverse problems. 
We show how the classical Birkhoff normal form can be recovered from semi-classical spectral 
invariants. In fact the full quantum Birkhoff normal form of the quantum Hamiltonian near a 
closed orbit, and infinitesimally with respect to the energy can be recovered. This generalizes 
recent results of Guillemin and Zelditch [|l^,|Q, |^ obtained in the high energy setting (a 
special case of semi-classical asymptotics). 

We will illustrate the results in a new setting to which they apply. Let P{h) be a semi-classical 
Schrodinger operator: 

(1.1) P{h) ^ -h'^A + V{x) - E , y G C°°(R"+^M) , 

with the principal symbol p{x, C) = + — E. We make the following assumptions: 

(1.2) 3e>0, p-i(he,e]) ^T*]R"+S dp ^ Q on p-^ {Q) . 

The first assumption guarantees that P{h) has a discrete spectrum near 0, and the second one, 
that the energy surface p~^(0) is smooth. 

The semi-classical inverse problem can be formulated as follows: 

What information about the energy surface, -I- V{x) = i?, can be recovered from 
the asymptotics of the spectrum of — /i^A + V{x) in a small fixed neighbourhood of 

E,a.sh-^Q ? 

This is very natural and can be considered as a mathematical formulation of the general prob- 
lem arising in spectroscopy. In concrete physical situations, the relation between the Hamiltonian 
and the spectrum is investigated using spectroscopic Hamiltonians which are essentially Hamil- 
tonians in Birkhoff normal forms - see for one of the first uses of Birkhoff normal forms in 
theoretical chemistry, and for recent applications using experimental spectral data. 

The energy surface, p^^(O) C T*M" has a natural bicharactistic foliation, coming from the 
Hamilton vector field. Hp — — dx^p d^^. Since it describes classical dynamics, 

its properties are a natural target of a recovery procedure based on quantum information, such 
as the spectrum. A naive intuition suggests that the quantum bound states (corresponding 
to eigenvalues) should correspond to the closed orbits of the classical dynamics. The simplest 
symplectic invariants associated to such closed orbits, 7, are the actions, 

— S^dx , 

the eigenvalues of the linear Poincare map, P-y, and the Maslov index, v-^ ^TL^. We recall that 
is the differential of the (non- linear) Poincare map, K-y, defined as follows: let C/ be a hypersurface 
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inp ^(0), transversal to 7, locally defined near some point on 7. Then [/ is a sympleetic manifold. 
The sympleetic map :[/—)■[/ is defined as the first return map: 

K^(m) = exp(T(m)iJp)(m) € U , eiip{tHp){m) ^ [/ , <t < T{m) , 

where, strictly speaking, the definition may require replacing [/ by a larger neighbourhood as the 
range of k^. 

The relation between the spectrum and closed orbits is provided by trace formulae: the 
Gutzwiller and Balian-Bloch trace formulae in physics and their mathematical precursor, the 
Selberg trace formula, and successors, the Duistermaat-Guillemin and the Guillemin-Melrose 
formulae (see [0 and references given there). Gutzwiller's formula states that if det(/ — P^) 7^ 0, 
for ^ kT^ e supp g, then for g,x ^ X = 1 near 0, g = near 0, 

(1.3) trg -pW(W)-E E M n.^' H^inAgW, a,.,^,^g{kT,). 

^ ^ 7 o#fc=-oo |aet(,i"^ 7j|2 ^.^g 

where is the period of 7, and the 7 summation is over all simple closed orbits of Hp on (0). 
The Maslov index, v^^k^ may in principle depend on the number of times of "going around" the 
orbit. The coefficients aj.^^u are distributions on R: C^(M) ^ g ^ flj,7,fc- 

It is now standard that under the non-degeneracy and simplicity assumptions, det(/ — P-y) 7^ 0, 
T~f — Tj' ^ 1 — l\ we can recover the actions, S^,, periods, T-y, Ma slov indices, z^^^fc, and 
|det(/ — P^)|, for all k. Under stronger non-degeneracy conditions (see (L4) and (2^) below), 
a theorem of Fried Q (see Proposition |3.2| below) then shows that the eigenvalues of the linear 
Poincare map, Py, can be recovered. 

The higher order coefficients aj^^,k contain further information and they form a family of semi- 
classical spectral invariants. It turns out that under a stronger version of the non-degeneracy 
hypothesis, we can recover from them the full infinitesimal information about the non-linear 
Poincare map n^. For simplicity we state it here in the case of elliptic orbits only, which is the 
semi-classical analogue of Guillemin's result - see Sect. I for the general statement. 

Thus suppose that the eigenvalues of dny^wo), wq — C/ H 7, are given by exp(±i0j), j — 
1, • • • ,n, 9j e (0,7r). We strengthen the non-degeneracy assumption to independence from 2tt 
over rationals: 

n 

(1.4) '^kj0jG2TrZ, kj €Z =^ kj = 0, j = !,■■■ ,n. 

The now classical normal form theorem of Birkhoff, Lewis and Sternberg says that there exist 
sympleetic coordinates (x,^) G M^", centered at wq, such that in a neighbourhood of wq, we have 

k^^ho + k\ k^(x,0 =0((a;'+e')°°), 

Ko{i,u;) = (^L,UJ- , t (tl, • • • , t„) , = (wi, • • • , w„) , 

n 

p{i) ^p{Q)-Y,e,L,+o{L^), 

see [0, Theorem 2.1]. The Taylor series of the twist map kq is called the Birkhoff normal form of 
K. It invariantly describes the infinitesimal properties of the fiow near 7. Since we chose 0j's in 
a unique way, Q < 9j < tt, the Birkhoff normal form is unique. 

Specialized to this case our general result Theorem ^ can be stated as 
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Theorem 1. Suppose that P{h) = -/i^A + V{x) - E is a semi- classical Schrddinger operator 
satisfying and that 7 is an elliptic closed orbit of Hp on p^^{0), with eigenvalues of P-^, 



e , satisfying (1.4). Then the leading term and the coefficients aj,j_k in (1.3), determine the 



Birkhoff normal form of 

As stated above, our motivation for the recovery of the classical Birkhoff normal forms came 
from an attempt to understand the results of GuilleminJ^ and Zelditch |l^,[^,|^, in the 
context of the recent trace formula of Sjostrand-Zworski [|17[, and of the (quantum) Birkhoff nor- 
mal form for semi-classical Fourier integral operator of lantchenko ||l^ and lantchenko-Sjostrand 
JI^ (see also (l6|). Guillemin and Zelditch considered the high energy (or C°° singularities) 
case, corresponding to h^Ag — 1 semi-classically. Guillemin's starting point was his earlier work 
with Prangoise Q , and an observation of Zelditch on the non-commutative residues of (classical) 
Fourier integral operator |^ . Zelditch's approach also used the non-commutative residue but was 
more concrete and computational. 

In addition to recovering the classical Birkhoff normal form, §, [0, |§, m were concerned 
with recovering the full quantum Birkhoff normal form of the operator near the closed trajectory. 
Without giving a definition of the Birkhoff normal form this result follows from 

Theorem 2. Suppose that the assumptions of Theorem [j a re satisfied for potentials V and V , 
with orbits 7 and 7, respectively. If, in the notation of ( |l.3| ), = T;y, and \ det(/ - 

Pi-')! = I det(/ - P|)|, v^^k = i>7,fc, for all k, then 



There exists a unitary h-Fourier integral operator F , 
microlocally defined in a neighbourhood of 7, 
such that {-h^A + V)F = F{-h^A + V), 
microlocally to infinite order along 7. 



The general result follows from Proposition 4.1 and Theorem ^ below, and is given in Corollary 



5.2. The definitions of /i- Fourier integral operators, and of agreement "to infinite order" are 
recalled in Sect. ^. 

Our approach is a straightforward application of (2.3) and ( 2.13| ) quoted below from our earlier 



papers, and of a result of Fried on the recovery of the eigenvalues of P-y (see Proposition 3.2 
below) . The trace formula shows that the traces of powers of the quantum monodromy operator 
of are spectral invariants. The quantum monodromy operator is a semi-classical operator 
which quan tizes the Poincare map for a given orbit 7, and once it is put into th e B irkhoff normal 
form using ( ^.13 ) the traces of its powers are easy to compute ~ see Proposition 3J . The recovery 



of the classical Birkhoff normal form is then very clear - see Theorem ^. 

One of the most striking applications of this approach to inverse problems is the result of 
Zelditch 111 on the recovery of bi-axial analytic planar domains from the spectrum of the (Dirich- 
let or Neumann) Laplacian. By applying our general methods we avoid any specific new work 
involving the boundary value problem - see Theorem ^ and Corollary |3.4[ . This and other inverse 
results are presented in Sect.||. 



What may seem like an excessive generality in the trace formulae of (see (2^) below) and 
Theorem ^ is motivated by their potential use for the recovery of classical dynamics of effective 
Hamiltonians. One of the most impressive inverse procedures is the application of the Onsager 
rules for determining Fermi surfaces - see for instance ^ for an introductory physics discussion. 
They can be interpreted in terms of a trace formula involving an effective Hamiltonian coming 
from the Peierls substitution - see Q . The question which will be investigated elsewhere is: can 
a more detailed information about the Fermi surface be obtained from higher order terms in the 
expansions of magnetic susceptibilities of metals. 
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Throughout the paper we will write 

oc 

to denote asymptotic expansions as e — > 0. When no confusion is likely to arise all equalities 
are meant modulo ^^(e"") where e is the relevant small parameter (mostly h, the mathematical 
"Planck constant"). We denote by N the set of natural numbers, 0, 1, 2, • • • . 

Acknowledgements. The third author would like to thank the National Science Foundation 
for partial support under grant DMS-9970614. He would also like to thank Bill Miller and Howard 
Taylor for references to the chemistry literature. 



2. Preliminaries 

In this section we will recall the general trace formula of jl7| and the classical and quantum 
Birkhoff normal forms of ||l^ . 

We consider X which is either a compact C°° manifold of dimension n + 1 or M"+^. 
We introduce the usual class of semi-classical symbols on X: 

S"''\T*X) = {aeC°°iT*X X (0,1]) : \d^d^a{x,^;h)\ < C'„,;3/i~"(0''"'^'} , 



and the corresponding class of pscudodifferential operators, 5'™''^(X), with the quantization and 
symbol maps: 

Op^J" : S"'-''{T*X) -^^';^^''{X) 

Gh ■■ 4'"'''(X) — > S'"^'=(T*X)/S'"-i''=~i(r*X), 

with both maps surjective, and the usual properties 

ah{AoB)=<Jh{A)<Jh{B), 

^ ^™-i.fc-i(x) ^''"■'^(X) ^ S"'^^{T*X)/S"'-^-''-^{T*X) , 

a short exact sequence, and 

ah o O^t : S'"'''{T*X) S"'-''{T* X)/ S"'-^^''-\T* X) , 

the natural projection map. 

Let P{z) E C°°{Iz','^'h'^{X)), I = (—a, a) C M, be a family of self-adjoint, principal type 
operators, such that Y.^ = {ni : (t{P(z)) = 0} C T*X is compact. We assume that 

a{d,P{z))<-C <Q, nearE,, \a{P{z))\ > 0(0" , 

\ ) 11 

for 1^1 > C when X is a compact manifold, and for |(a;, 01 > C if X = R""*" . 

We also assume that for z near 0, the Hamilton vector field, Hp^^^^ p{z) — a{P{z)), has a simple 
closed orbit 7(2;) C with period T{z), and that 7(2) has a neighbourhood fl such that 

(2.2) 

m Q and exp tHp(^z-^ (m) = m , p(z, m) = 0, Q < \t\ < T {z)N + e , z Cz I , m G 7(2) , 

where T{z) is the period of 7(2), assumed to depend smoothly on z. We also write p{z,m) for 
a{P{z)){m). Let A G 5'°'°(X) be a microlocal cut-off to a sufficiently small neighbourhood of 
7(0). 
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Then if P(z) is an almost analytic extension of P{z), z e M, x e C^{I), x G C^{C), its almost 
analytic extension, / e C°°(R), and supp / C i-N{Cp - e) + C, N{Cp - e) - C) \ {0}, we have, 

itr J f{z/h)d, [x{z) d,P{z) P(z)-i] AC{dz) = 
(2.3) N-i 



tr / f{z/h)M{z,hf^M{z,h)x{z)dz + 0{h^), 

where M{z^ h) is the quantum monodromy operator. The constant Cp > 0, in the condition on 
/ depends on p(z). 

The quantum monodromy operator, M{z, h) is defined as follows. For a point on 7, mo G 7, we 
can define the microlocal kernel of P{z) at mo, to be the set of families m(/i), such that u{h) are 
microlocally defined near mo and P{z)u{h) = 0{h°°) near mo- We denote it by ker„(j(2-)(P(z)). 
Any solution can be continued microlocally along 7(2) and we denote the corresponding forward 
and backward continuation by I±{z) - see §4] for precise definitions. Now, let mi ^ mo be 
another point on 7(2). We then define 

I^(z)Miz) — I<(z) , near mi, 

(2 4) 

M{z) : ker„„(^) (P(z)) — > ker„,„(^) (P(z)) . 

The operator P{z) is assumed to be self-adjoint with respect to some inner product (•,•), 
and we define the quantum flux norm on ker„jfj(2) (P(z)) as follows: let x be a microlocal cut-off 
function supported near 7 and equal to one near the part of 7 between mo and mi. We denote 
by [P(z), x]+ the part of the commutator supported in mo, and put 

(w, u)qf =^ {[{h/i)P{z),x]+u,v) , u,v & ker„^(^)(P(z)) . 



It is easy to check that this norm is independent of the choice of x ~ see |17, Lemma 4.4]. This 
independence leads to the unitarity of Ai{z): 

{M{z)u,M{z)u)q-p = (m,u)qf , u G ker„„(^)(P(z)) . 

For practical reasons we identify \ieYmo{z){P{z)) with I?'(R"), microlocally near (0, 0), and choose 
the identification so that the corresponding monodromy map is unitary (microlocally near (0, 0) 
where (0, 0) corresponds to the closed orbit intersecting a transversal identified with T*R"). This 
gives 

M(z,/i) : 2?'(R") — > V'{W), 

microlocally defined near (0, 0) and 

M(z, /i) is a semi-classical Fourier integral operator which quantizes the Poincare map 
of 7(2;). This, and (pT^), ( |2.6| ) are the only properties of M(z,h) which will be used 
in this paper. 

Basic ideas become clear when one considers the simplest example P(z) ~ hD^ — z, on S^. In 
that case M{z, h) = exp(27riz//i) - see [0 §2] for a careful presentation. 

We will now assume that the symbol of P(z) has an asymptotic expansion in powers of ft., 
near Eg: 

(2.5) 

00 

{P(z)){x,S,\h) —'Y^y{pj{z)){x,C)^ i^iO G a precompact neighbourhood of So . 

3=0 
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Under this assumption wc have a simple extension of [ pT| Proposition 7.5] 



(2.6) trM(z, hf-^hD,M{z, h) = e'''^'-'^ ^ + e''''^'^ ^ aj^kh' , 

j=i 

where I{z) is the classical action, the Maslov index of the fcth iterate of 7(2) (it is independent 
of z for z near 0). 

We start by recalling the general statement of the classical Birkhoff normal form - see |p^ , 
Theorem 1.3, Proposition 4.3] and references there. 

Let (VF, w) be a symplectic manifold of dimension 2n, and let k : W — > W , h{wq) — wq , 
be a symplectic transformation, k*u> = u, fixing wq G W. In (1^), the twist map, kq, is the 
Birkhoff normal form of k, 

(2.7) Ko = exp_ffp, 
where the Hamiltonian p is in the Birkhoff normal form. 



The general case presented in [g^ is formulated using (2/7): suppose that the eigenvalues of 
dn{wo) are given by 

X] ,Kj^ , j = !,■■■ ,n 

Xj — exp , fij = aj + 1(3 j , > , < fij ^ 2ttZ , I < j < rich , 

(2.8) Aj+„^i, = exp/ij+„^j, , Aij+nch = 7*7, 1 < j < "-ch , 

Xj = exp fij , /ij > , 2nch + 1 < J < 2nch + "-rh , 

Aj — exp fij , /ij e i(R \ 27rZ) , n^h + 2nch + ^ l£ j l£ n . 

Here ne + ?irh + 2nch = and e stands for elliptic, rh for real hyperbolic, and ch, for complex 
hyperbolic. 

We replace the condition in the purely elliptic case (|l.4D by the more general condition 



(2.9) kjfij e 27riZ , kj e Z =^ = . 

i=i 

Then the generalization of ( |l.5D says that there exist symplectic coordinates (.t, ^) e M" centered 
at (0,0) 

(2.10) 

K = exp Hp + k\ k\x,0^O{{x^ +^^)°°), 
p = Po + r, r{x,0=O{{x'+e)'), 

"ch 

Poix,0 = +X2j£,23) - Pj{x2j~1^2j " X2j^2j-l) 

2nch+"rh " 

r(x,^) = i?(i), i?(0= E 

a|>2 

(a;2j-i6i-i + 2;2jC2j + iix2j^i^2j - a;2jC2j-i))/2 1 < j < "ch 

Lj> j' = j - rich , rich + l < j < 2nch 

a^jCi 2nch + 1 < .7 < 2nch + ^rh 

(2;? + ^?)/(2*) 2neh + Tlrh + l<3 <n. 
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Observe that pq — (t, /i). By the Birkhoff normal form of n we mean po and the asymptotic 
expansion of r. 

We now recall the quantum Birkhoff normal forms [|ll] and ||l2|, Theorem 5.1, (4.11)]. We start 
with some general facts about /i- Fourier integral operators. Let W above be a neighbourhood of 
Wo S T*X, where X is a C°° manifold. To k : W ^ W we can associate a class of operators, 
microlocally defined near wq (see |l^, §3] for the precise definition of this notion): 

the semi-classical Fourier integral operators quantizing k. One way to define them is by consid- 
ering oscillatory integrals: we identify wq with (0,0) S T*R" so that 

U e riC{X,X,K') Uuix)^!!-"'-"^ J J e"l'<-^'y-'^'^/''a{x,y,9;h)uiy)dyd9, 

(2-11) 0eC°°(M"xM"xR^), graph(K) = 7r(C0) , 

C4, = {{x,y,e) : 4>'g{x,y,0) ^ 0} 3 {x,y,0) {x,(j>'^;y,(j)'y) , 

where is assumed to be a non-degenerate phase function in the sense of Hormander (see |jl^, 
Dcf. 21.2.5] for a related discussion), that is (p is smooth and real, is a smooth manifold of 
dimension 2n, and tt has an injective differential. The am plit ude a is assumed to be supported 
near (0,0,0), and to be a classical symbol in the sense of (2.5). 

Before stating the main result of ]|l^ we have to introduce the notion of equivalence of families 
of operators used in that paper H §2, §5]: let U{z) e l"{X,X, 4) and U{z) e l"{X,X, k'J be 
two families of operators. Then 

(2.12) U=U 

if ^0(^0) = Ko(u'o) = Wq, the two families of transformations, Kz and k^, agree to infinite order 
at Wq and z — 0, and the terms in the asymptotic expansions (in powers of h) of (jjz, 4>zi o-z, dz 
(with amplitudes as in ( 2.11 )) agree to infinite order at (0,0,0) and z = 0. 

Suppose that U{z) quantizes Kz and kq satisfies the assumptions of ( ^.10[ ). Also assume that 
U{z) is elliptic. Then 

U{z) = V{z)-^e'P^''''^/''V{z), 
V{z) e l'^{X,X,Cz) is a family of elliptic Fourier integral operators near ((0,0), wq) 
P{z, h) = pq{z) + hpi{z) + h?p2{z) + • • • , 
Pj{z) =pj{z;Ii,-- ■ ,Ij) , Ij = Lj{x,hD.j;) , 



(2.13) 



When U{z) is microlocally unitary near wq, P{z, h), V{z) can be chosen to be microlocally self- 
adjoint and unitary, respectively, at ((0,0), wq). We will call the expansion of P{z,h) at z = 0, 
L — 0, the infinitesimal Quantum Birkhoff Normal Form of U{z) at z = 0. 

We recall that a Fourier integral operator is elliptic if it is associated to a canonical transfor- 
mation, that is C above is the graph of a canonical transformation, ((0,0), wq) G C, and that its 
symbol is elliptic, that is the leading term in the asymptotic expansion of the amplitude in (2.11) 
does not vanish. 
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3. Recovering Birkhoff normal forms from traces 

From Lemma 3.2, Proposition 7.3] we know that if U{z) G I^{X, X, kIz)'), is a family 
of /i-Fourier integral operators, and that the canonical relations n{z)^ , k ^ 0, have only one 
non-degenerate fixed point then 

(3.1) tr Uiz)"" = e'kS{.)/h ^ ^^^(^^/^ ^ . . . ^ ^^^^^^/^i ^ 0{h'+^)) . 
In this section we will prove the following 

Theorem 3. Let U{z) e I^^(X , X , k(z)') be a family of h- Fourier integral operator associated to 
smooth locally defined canonical transformations, k{z) : T*X — s- T*X. Suppose that kq has a 
unique fixed point, wq, and that the eigenvalues of dKoiwo) given in (2.8) satisfy (2.9). Then 
the Taylor series of the coefficients aik{z) in (3.1) determine the infinitesimal quantum Birkhoff 
normal form of U (z), P{z, h) in (2.13). In particular, the Birkhoff normal form of kq at z = 0, 
p{z), in ( 2.1C ), is determined. 

Our proof of Theorem ^ is based on 
Proposition 3.1. Suppose that 

(3.2) G{i;h)^{L,ii)+Fii,h), FiL,h)^FoiL) + hFiii) + -- - , e^' ^ I , Fo{i) ^ 0(1,^) , 

where we consider F(l, h) as a formal power series in h, and each Fj{i), a formal power series 
in L, and where we supressed the dependence on fj, in higher order terms. 
Iflj{h) = Lj{x,hD^), then 

(o o^ ^ -iG{I{h)-h)/h _ -iF{ihd^;h)/h I A ^ 

" Ul2sinh(^,/2) 



Proof. When G{l; h) = this follows from computations in each of the cases given in ( 2.1C ): 

where we can separate the trace into separate blocks: 
• elliptic blocks: 



oo 



2 sinh 



where the formal computation (/i G iM+) can be justified by condsidering the limit e — > 0, 
with /Lt replaced by /i + e. 
real hyperbolic blocks: 

Formally, 

and this can be justified by approximating 6q by a sequence of smooth compactly supported 
functions. 
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• complex hyperbolic blocks: here we have to work in two dhiiensions, /ii = a + i/?, /i2 = /ii = 
a — i(3, and the corresponding ii and L2 (see the last two hnes in ( [2.10[ )) are ti = z • 
L2 = w • ui, where z = {xi — 1x2)/ V2, w — z and C = (^1 + «^2)/ \/2, = (Ci ~ *'^2)/ \/2, the 
corresponding dual coordinates, uo — C,- With Ij = lJ{x, D^), j = 1, 2, we have 

The last operator can be written as a contour integral, 

with the integration contour given hy w — z , lu = 

Formal stationary phase argument applied when taking the trace gives 

n— 

2sinh^ ■ 

j=i ^ 

Once we understand the contribution of the leading term the following formal argument es- 
sentially gives ( t3.3[) . Since I{h) is unitarily equivalent to hl{l) we see that 

tr ^-^G(I(h);h)/h ^^^^ ^-lG{hI(l);h)/h _ 

On the other hand, /(I) = 9^(7(1), /i), and hence 

^-iG{hI(l);h)/h ^ ^-iF(hI{iyji)/h^-i{I(l)4L) ^ ^-iF{ihd^,;h)/h^-i{I{l),f,) ^ 



By moving differentiation outside of the trace we obtain the expansion (B.3). □ 

We also need a result of Fried ||^ used in the same context in : 
Proposition 3.2. // the assumption (|2.9|) holds then the set 



]^sinh(fc^j/2) , ken, 
i=i 

determines exp fij 's and gq uniquely. 



□ 



Proposition 3.2 and an argument based on Kronecker's Theorem give the following inverse 
result crucial in the proof of Theorem || (see 0, §8], |20, §6]) for different original approaches): 



Lemma 3.3. Suppose that fij satisfy (2.9). Then for any polynomial, p{^) = '^\a\<Ai '^a^" 7 
coefficients Oq, can be recovered from the asymptotics of 

n 

^(-^'^^-)n 2.inh(fc,,-/2) ' 



Proof. Using the argument of the proof of Proposition p . l| we can rewrite (3.4) with cq = (1, 1): 

^3 g-j j=l V f*J/ / mGN" 

= J2 p{{m + eo/2)/i)e~^"'+'°/^'''>'K 

meN" 

To justify this expression we will consider it as a distribution obtained as a boundary value of a 
holomorphic function. We introduce new complex variables z S fii = £^(0, 1)"° and w G il,2 = 
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{D{0, 1)\0)"'^'' in a polydisc and a punctured polydisc respectively: at the boundary, zj — e'^^^^, 
and Wj — e*'^J+"'=+"rh/^. We then write 

m = (m^^h, Kh, rrirh, m^) G N" , cq = {ef", ef", ej,'', eg) , 

^ = (^ch, A^ch, Mrh, Me) , ^J.ch = a + ip , ^i-h e R+ , ^J■c e , 

(m + eo/2, /i) (n2 + ef", a) + (m + e^'V2, /i^h) + i{h + e72, ^c/i) + iih, P) , 
Here m € N" is arbitrary, and 
with the constraint that 

2 

which imphes that IZ2I < ?^2, component wise. We then put 

9(?i,0 =p((m + eo/2)/i), z = e"^" G C"" , w = e^^'^/^ G C"<=\ 



(component-wise in the natural sense). In this notation (3.5) becomes 

(3.6) e-^«"i+^o"/2.A'.0 + {"2+eS'72,a»^^^(^fe^y,fe)^ 

ngN"rh + "ch 

where 

F„(z,H= 5] g(n,0z2'^+^°«)''% 

gN"ch 

(z, w) G fil X C C"'' X C"'^''. Here we note that for any fixed n, g(n, Z) is non-zero for finitely 
many Za's and hence Fn{z,w) is holomorphic in the punctured polydisc. Our task to recover 
the coefficients q{n, I) from asymptotics as fc — > 0. That amounts to recovering the holomorphic 



functions F„. We observe that (3^) shows that the distributional boundary values F„'s are 
smooth away from zj = 1. 

An application of Kronecker's theorem]^ shows that for any 

(a;,y)G(9i?(0,l))"<=x(9i?(0,l))"-, 

there exist a sequence ^ 00 such that (z'^'', u'*^'^) (a;,y). 
From ( |2.9| ) we deduce that 

kj Re /ij = kj — . 

1=1 

Henc e by choosing sequences of fc ^ 00, we can determine the coefficients, Fn{x,y) in the first 



sum (3.6) (the terms have different rates of exponential decay). From that we determine the 



holomorphic functions Fn, and the coefficients, q{n,l). □ 

Proof of Theorem^. Applying ||l|. Theorem 3.2] and Theorem 4.4] we see that, up to 0(/i°°) 
errors, we can write 

tr U{zf = tr e-k^G(i(h),z;h)/h ^ 



-"-which says that if a G Z" satisfies {a, k) G 27rZ, k € Z", fc = 0, then f or a ny y g R" and e > 0, there 
exists Af 6 N and K eZ", such that \y ~ Na- 2nK\ < e. Its use is justified by (p!s|). 
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where G is as in Proposition Hence, 

tr U{z)^ _ ^-ikF{ihk-^d^,,z:h)/h ^-i{I(l),k^i) 



Recovering the BirkhofF normal form of k means recovering the Taylor series of {i, i-i{z)) + 
Fq{l, z) from the trace. We will actually recover the full infinitesimal quantum Birkhoff normal 
form, ((., + F{l, z; h). 

We write 



F{hy,z;h) 



oo i+1 



3=0 



1=0 \a\=j+l-l 

oo 

/j(^.y) -E•^J■fe(2/)^''■ 



A;=0 



The polynomials fjk{y) clearly determine the formal series F{i,z;h). 
Since 

oo 

e^p{—ikF{hik^^df^,z,;h)/h) — e:x.p{—ik'y^h^fj{z,ik^^dfj_)) = 
1 - i ^ h^kf,{z, ik-^d^) -ttAY. h'kfA^^ 



-ikfo{z) 



\ 



i=i 



we see from (|3.3| ), that the coefficient of in the expansion of the trace is given by 
e-^'^^™ ( (1 - * E ^"^/o™ ~ ^ f £ ^"fc/o™) + • • • ) X 

\ m— 1 ' \m— 1 / / 

/ \ 2 

l-^Y^zOk{^k-^^^,^Ji^=\0))-- [Y.z^k{^k-^^^,^i^^\0)) +• 
i=i ■ \j=i J 



tr e 



-j(7(l),fcM>| 



where we expanded into its Taylor series and then used the same argument as in the proof 
of Proposition 3.1. We apply Proposition 3.2 to the constant term to recover expi/oo and /i(0). 
For j > 0, the coefficient of z^ in the expression above is given by 



-ike-^'f^ih, + {k-hd^,fi^^\0)) + H,^k{ik-^d^)) JJ 



^ijL 2sinh(fc^j/2) I Im=p(0)' 



where is constructed from /o/'s and /.t^^-'(O) with l,p < m. Hence Lemma 3.5 shows that 
/oj's and y^^^'(O) can be successively recovered. 

We now move to higher powers of h, with the procedure being essentially the same: we see 
from (^^) that the coefficient of ft,™ in the expansion of the trace is given by 



e-'''M^\-ikUiik-'df,, z) + GmM^k-'d^n ^)) n 



2sinh(fc^j/2) I ' 

where Gm.k is constructed from /j's with I < m, and we already know the z-expansions of fo{z) 
and fi{z). Replacing by its expansion in z, reduces the problem to recovering the polynomials 
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fmiiy) from 



n 1 



keN, 



and that follows again from Lemma 3.3. This shows that we recovered the full expansion (in l, 
z, and h) of G(i, z; h), that is the infinitesimal quantum Birkhoff normal form at z = 0. D 



4. The monodromy operator determines P near 7. 

The trace formula (2^) and the results of Sect.^ will be used to recover the infinitesimal 
quantum Birkhoff normal form of the quantum monodromy operator M{z) at z = 0. In this 
section we will discuss the recovery of the full Hamiltonian P{z), infinitesima lly at z = 0, from 
M(z). Since our information comes from the left hand side of the trace formula ( |2.6| ), this recovery 
can only be possible up to conjugation by families of elliptic /i-Fourier integral operators, and up 



to multiplication by families of elliptic /i-pseudo-differential operators. The assumptions (2.1 
the implicit function theorem, and the usual symbolic iteration, imply that 

P{z)^Aizr{P-z)A{z), 

with A{z) e *°'''/^(X) ehiptic near 7(0), and P e *°'°(X) is self-adjoint. Replacing P(z) by 
P — z in ( |2.6| ) changes the trace by 0{h°°). Hence we would like to recover P infinitesimally near 
7, on the energy surface p~^(0), p = cr(P). 

Thus in this section, let P be a self-adjoint ft,-pseudodifferential operator with leading symbol p. 
Let 7 be a simple closed i?p-trajectory of period Tq > 0. Recall from |l^, §4] that the quantum 
time, Q, is an ft,-pseudodifferential operator microlocally defined in a suitable neighbourhood 
W+ = neigh (toq) of some fixed point toq G 7i such that 

'-[P,Q]=Id. 

After conjugation by a unitary Fourier integral operator, we may assume that 

(4.1) P = hD,^, Q^xn. 

Let K{z) : I?'(R") V^X) be a microlocally defined Fourier integral operator with 

(4.2) (P-z)X(z) = 0, '-K[zY[P,x\w^K{z)^Id, 

where x G C°°(W+) is equal to near the incoming part of 7 and equal to 1 near the outgoing 
part, and we let x s-lso denote a corresponding /i-quantization. In the case of (4.1), we can take 

(4.3) K{z)u{x) = e^'^-Z^uix'). 



It follows from the second part of (4.2), that 

K{zrf'{Q)K{z) = Id, 

for / g C°°(R), = for large negative t and — 1 for large positive t. Indeed, e^'^/^'Qe"**^/'' 
Q + t and 

J[P,/(Q)] = (50t=o(e^*'^/V(Q)e-'*^/") = (9t)t=o/(e^*''/'Qe-'*^/'') 
^{dt)t=ofiQ + t) = r{Q). 
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Notice that /' e C^^'CR), J f'{q)dq = 1. More generally, 

(4.4) K{zrg{Q)K{z) = j g{q)dq ■ Id 

for g G C^(R). As in Proposition 4.3] we can choose K{z) so that 

(4.5) K{z) = e*(^-"')'3/''X(u)) , 

1) hold (note that for P{z) = P — z, U{z,w) = exp(i(z — w)(5), 



since this is fulfilled when ( |4.l[ 
in (4.4)]). 
Consider 

^{z)"^' ^K{z)K{zY. 

Notice that this definition is independent of the choice of K{z) satisfying ( |4.2| ). 
We have the properties: 



7r(z) = 7r(z)*, 
P7r(z) = 7r(z)P = Z7r(z), 

T:{z)dz = Id , 



(4.6) 
(4.7) 

(4.8) 

where for the understanding of ( |4.8| ), we remark that if S C^(R) is independent of /i, then 
/ 4>{z)T:[z)dz is a well defined /i-pseudodifferential operator of order 0, and that we can always 
choose (j) to be equal to 1 in some interval corresponding to the region of phase space where we 
wish to establish (4.1 



The relations (4.6), (4.7) are obvious, and so is (4.8) if we work in the representation (4.1), 



We also give a direct proof of (4.8) which docs not involve conjugation to the model: Let 
g € C(5"(R) and consider 

K{zyg{Q)K{w), 



where we will assume (4.5). Then we get 

(4.9) K{zyg{Q)K{w) = K* {z)g{Q)e-'^'-^'^Q"^ K{z) = {Thg){z ^ w) ■ Id, 

where J^hg{z) — J e^^'''^/'^g{q)dq is the /i-Fourier transform. Here the second equality follows 
from {iA). 

Consider the microlocally defined /i-pseudodifferential operator: 



A 



n{z)dz. 



A is an elliptic pseudodifferential operator of order 0, and we claim that 

(4.10) A^ = A. 
Since the statement is microfocal, it suffices to show that 

(4.11) Ag{Q)A = AgiQ), 



g e Cff (and then choose g{Q) = 1 in the region of interest). To get (4.11), we first compute: 

1 



'^iz)g{Q)'K{w) 



;K{z)K{z)*g{Q)K{w)K{w)* 



^-L_K{z)K{zr{Tug){z - w)e^^^--)Q'\ 
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where we used (|4.9|), ([4. 41). Integrating with respect to w, we get 

7r{z)g{Q)A^7T{z)g{Q), 



which gives (4.10) after integration in z. 

Rewriting ( 4.10| ) as A{A — 1) = and using that A is elhptic, we get 



Remark. Using (J^^ (or a direct argument), one can check that 

where the integral to the right is microlocally well-defined by restricting the integration to a 
suitable finite interval. 

Let P, P be two ft,-pseudodifferential operators with simple closed trajectories 

7Cp-i(0), 7Cr'(0) 

with the same period Tq. Let ttiq G 7, ffiQ e 7 and assume that the corresponding monodromy 
operators 

M (z) : kermo {P - z) — > ker,„„ (P - z) , 
M {z) : keirho (P - z) — > ker^^ [P - z) , 

are conjugated: 

(4.12) Uiz)-'^M{z)U{z)=M{z), 

where U{z) : kermp(P — z) — > 'keTmoiP — z) is a unitary Fourier integral operator depending 
smoothly on z. 

Identifying P, P with /il?x„ , we see that 

F =^ f U{z)Ti(z)dz 



is an elliptic unitary Fourier integral operator with 

We can identify P, P with hD^^ in such a way that F becomes the identity operator. With 



this identification, we take K{z) — K{z) of the form (L3). Without using the identification, this 
means that K{z) = FK{z) and 

(4.13) ^ = i / K{z)K{zrdz. 
Also, 

(4.14) M{z) = M{z). 

Let Kf{z), Kf{z); Ki,{z), Ki,{z) be the forward and backward extensions of K; K as in |lj, §5]. 

If u is microlocally concentrated near m = exp(ii?p)(mo), — e < i < Pq ^ 2e, where Tq is 
the common period of 7,7, we define P/it, microlocally concentrated to a neighbourhood of 
exp(iiJp(mo)), by 

(4.15) Ffu = e-'^^'^Fe'^P'^u. 
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Since F intertwines P and P, this definition is invariant under small variations of t and we see 
that Ff is a unitary Fourier integral operator with canonical transformation, given by 

KFjijn) — exp(iiJp) o Kp o exp{—tHp){m) , 

for m in a neighbourhood of exp(t_ffp)(mo), and where Kp ■ neigh (mo) neigh (mo) is the 
canonical transformation associated to F. 
We have 

Pf^^^l Kf{z)Kf{zydz. 



In fact, by ( [4.131 ), ( |4.15| ) 

^ 2TTh 
__ 1 

~ 2TTh 



e-~'''/^Kf{z)Kf{z)*e'^''''dz = 



1 

2TTh 



Kfiz)Kf{zrdz. 



Similarly, if u is microlocally concentrated near m = exp(<i7p)(mo), — Tq + 2e < t < e, we 
define FbU to be the right hand side of (4.15|) (now with t in a different interval) and get 



Fbu^-^ I Kb{z)Kb{z) 
Znn 



dz. 



If m belongs to the overlap region, where both the forward and the backward extensions of F 
are defined, then for u concentrated to a neighbourhood of m, we have FfU — F\,u: 



1 



1 



FfU = — / Kf{z)Kf{z)*udz = I Kb{z)M{z)M{z)*Kb{z)*udz ^ Ftu, 



2nh 



since M{z)M{z)* = M{z)M{z)-^ = 1, in view of ( |4.14| ). 

This means that F extends to a well-defined operator in a neighbourhood of 7, and we get 

Proposition 4.1. Let P, P be two h-pseudodifferential operators with closed simple trajectories 
7 C p^^(O), 7 C p^^(O) having the same period. Let mo G 7, too S 1, o,nd assume that the 
monodromy operators 

M{z) : ker„, (P - z) — > ker„„ (P - z) 
M(z) : ker^o {P - z) — > ker„-i, (P - z) 



are conjugate as in (4.12), via an h-Fourier integral operator U{z) depending smoothly on z. 
Then there exists a unitary Fourier integral operator F , microlocally defined in a neighbourhood 
of J with an associated canonical transformation Kp : neigh (7) — > neigh (7), such that PF = FP . 

□ 

Using Theorem ^ we will recover t he in finitesimal quantum Birkhoff normal form of A^(z) at 
z — Q. From the proof of Proposition 4T it is clear that the conclusions remain valid to infinite 
order at 7 and 2: = 0, if that assumption is made on the conjugation. Hence we obtain 

Proposition 4.2. Suppose that in Proposition 4-1 we only assume that 

M{z)=U{z) M{z)U{z)-^ 



to infinite order at m — toq, z — Q in the sense of (2.12). Then there exists a unitary Fourier 
integral operator F , microlocally defined in a neighbourhood of 7 with an associated canonical 
transformation k p : neigh (7) — > neigh (7) , such that 



PF — FP , to infinite order along 7 . 



□ 
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5. Applications to inverse problems 



We start with a simple lemma motivated by ( p.3[ ) and (2.6): 
Lemma 5.1. Let u{z, h) be a semiclassical family of functions defined near z — 0, 

oc 

u(z,/i) = e*-^(^)/''^aJ(z)/^^ /, G C°° , ao(0)7^0, lm/(z) = 0. 

3=0 

Suppose that V C C^' (M) is a subspace from which we can pick an element having any prescribed 
finite Taylor polynomial at I'{0). If for any g Cz V we know bj{g) where 

/oo 
g{zlh)x{z)u{z, h)dz = e^^^")/" ^ bk{g)h'' , x G C^W , X=^ near 0, 

k=0 

then the distributions bk determine the /''^(O), and a^p (0), for all I and j. 

Proof. We first observe that if aj{z) = 0{z°°) then they contribute 0{h°°). Hence in considering 
the expansion of J{g, u) we can replace / and a^'s by their formal power expansions: 

C30 oo 

I{z) = lo+hz + Y, h+iz'"^^ , aj{z) = '^okz'" , 

fe=l k=0 

and our task is the recovery of /^'s and ajk^s. Making a change of variables C, = z/h we write 
J((7, u) as the formal integral 



where we can neglect x(/iC)i ^ x[z) = 1 near z — Q. Just as in the proof of Theorem^ we now 
see that the coefficient of Kp^^ in the expansion is given by 

e^"/"/! j g{Cy''^ ( a;feC'+*aoo/p+iC^+'+i?™K/) I , 

•' \l+k=p ) 

where Hm{a, I) depends only on aik with I + k < p, and Ii with I < p + 1. Since we assumed 
that ooo = ao(0) 7^ 0; ^^nd that we can find g £ V with any prescribed finite Taylor polynomial 
at Ji, the coeffients of / and a can be recovered from bj{g)'s, g G V. □ 

Using Lemma 5.1, and Theorem 2, Proposition 7.5] (see (|U) (U) above) we obtain the 
following consequence of Theorem ^: 

Theorem 4. Suppose that P{z) satisfies the assumptions needed for (2.3) and is classical in 
the sense of ( ^.5| ). Suppose that in addition the eigenvalues of the linear Poincare map of j{0) 
satisfy (2.9). Then the coefficients in the asymptotic expansion of 

_ 00 
-tr f{z/h)d4xiz)d,P{z)P{z)-']ACidz)= Y e'*^^/''^a,■fc(/)/^^ 
^ fcez\{o} j=o 

determine the infinitesimal quantum Birkhoff normal form of the quantum monodromy operator 
M{z) at z = 0. In particular the Birkhoff normal form of the Poincare map 0/7(0) is determined. 
Here the notation is the same as in (2.3), and by the coefficients we mean the distributions, 
/^a,(/), /GC-(K\{0}). 
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Proof. We first see that the non-degeneracy of j{z) for z close to shows that (2.6) holds. We 
apply and integrate by parts: 

1 r d 

J2 tr f{z/h)Miz,h)'' — M{z,h)x{z)dz^ 

h~'— f[z/h){k + l)-Hi{z,hf+'x{z) + 0{hn. 

Zm ^ — ' TO 

where x = 1 near and hence the contribution of x' is negligible: f{z/h) — 0{h°°) for |z| > e > 0. 



The assumptions of Lemma 5.1 with g ^ f are satisfied and hence we obtain expansions of 
trM(z, ft,)*^, for any k. Now we can apply Theorem || with U{z) = M{z,h) and the proof is 
completed. □ 



When we combine this result with Proposition 4.2 we obtain infinitesimal information about 
P[z) at 7: 

Corollary 5.2. Suppose that the assumptions of Theorem ^ are satisfied and P{z) — A{z)*{P ~ 
z)A{z), where A{z) is elliptic near 7. Then, up to conjugation by elliptic h-Fourier integral 



operators, in the sense of Proposition 1. 1, P is determined to infinite order at 7. 
When specialized to the classical, or high energy, case, we obtain 

Corollary 5.3. Let M he a compact C°° manifold and let P he a positive, self-adjoint elliptic 

1 

classical pseudo- differential operator acting on C°°{M, Suppose that 7 C p^^(l) is a simple 

non-degenerate orbit of the classical flow of the principal symbol of P, and P^ is its linear Poincare 
map. If ( |2.9D holds for the eigenvalues of P^ ( |2.8| ), and no other orbit inp~^{l) has the period 
which is a multiple of the period of 7, then the spectrum of P determines the Birkhoff normal 
form of the Poincare first return map. 

Proof. Suppose that the order of P is m. We define a semi-classical operator, P{h), by putting, 
P{h) = h^P{x, Dx) — 1 , so that the (semi-classical) principal symbol of P{h) is the same as the 
(classical) principal symbol of P. The non-degeneracy assumption for 7 C {p = 1} implies that 
the assumptions of Theorem ^ are satisfied, with P{z) = P{h) — z. The simplicity of 7 then 
shows that for /'s with supp / compact, and sufficiently close to ZT^, where T^ is the (primitive) 
period of 7 on {p = 1}, we have 

(5.1) tr f{P{h)/h)x{P{h))A = tr f{P{h)/h)x{P{h)) + Oih^ . 

The left hand side is a (semi-classical) spectral invariant and hence the result follows from The- 
orem ^. □ 

We leave it to the reader to formulate the corresponding corollary in the more general semi- 
classical case. 

In the case of manifolds with boundaries we consider, for the sake of simplicity, the case of 
second order differential operators only: 

Theorem 5. Suppose that {M,g) is a Riemannian manifold with boundary, and that Ag is its 
Dirichlet or Neumann Laplacian. 

Let ^ be a simple non- degenerate orbit of the classical billiard flow of g. If all the intersections 



of "f with dM are transversal, (2.9) holds for eigenvalues of the Poincare map of j (2.8), and 
the lengths of other closed orbits of the generalized billiard flow do not accumulate at a multiple 
of the period of ^ , then the spectrum of Ag determines the Birkhoff normal form of the Poincare 
first return map ofj. 
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Outline of the proof: First we need to indicate why the results of |T7[ can be applied in this 
situation. For that we will check that local solutions can be microlocally defined near a trajectory 
7. Let p be the symbol of — Ag and q the defining function of dM. The transversality condition 
on 7 means that Hqp 7^ at the boundary. We first recall (see for instance Sect. 24. 2]) that 
we can find coordinates in which, near a point in dM, q = xi, and p{x,£,) = Ci ~ ''(^^jOj 
^' = (^2, • ■ • jCm+i)? = dim M. Solving the boundary problem microlocally is equivalent to a 



hyperbolic problem and hence we have a solution. For localization needed in (2.3) we can use, 
near the boundary, an operator A = A{x, hD^'). |^ 

A more subtle argument is needed for showing that ( ^.l|) holds, as that involves results on 
propagation of singularities for boundary value problems - see [|o[ Sect. 24. 5] and references 
given there. These results provide the foundation for the trace formula of Guillemin and Melrose 
(see also |l^), which among other things shows that singsupp (tr cos(t^— Ag)) C £, where 
£ is the set of length of closed orbits of the billiard flow. It also follows from |^ and 
that if A{x, {Dt)~'^Dx) is an operator which provides a localization to a neighbourhood of the 
bicharacteristic strip of 7, then 

(5.2) trcos(tv/-Ag) - tr cos{ty/~ A g)A{x, {Dt)^^D^) e C°°(R) , for t near fcT^, 



that is, the singularities of the trace of cos{tyj'—Ag) at kT^ come from a small neighbourhood 
of 7 (depending on k). As stated above we only need tangential pseudo-differential operators to 
localize near 7. 



A straightforward translation to the semi-classical setting (see for instance |23|) gives (5.1). 
In fact, by standard oscillatory testing, A{x, {Dt)~^Dx)e^*^^ = e'*/''A(x, hD^; h) , where A has 
the semi-classical principal symbol given by A{x,£^). This shows that for / e C^(R) even, 

tr f{{h^/^g - l)/h)A{x,hD.,,h) ^ ^ J tr {A{x,{Dt)-^D,)co^t^f^^ f{t)e''"^dt , 



and thus (5.2) follows from (5.1). □ 



By combining Theorem |J, or rather Corollary 5_^, with Proposition 4^ or with the geometric 
arguments of §2] we obtain the Weinstein conjecture first proved by Guillemin in the case 
of elliptic trajectories, and by Zelditch [l9| ,[pot in general. It involves the concept of the Birkhoff 
normal form for the Hamiltonian - [fzL Theorem 1.1]. As shown by Frangoise and Guillemin 



and Guillemin |Q §2], and as can be also deduced from Sect.^ above, that normal form is 
essentially equivalent to the normal form of the Poincare map. 



Corollary 5.4. Let (Af , g) he a compact Riemannian manifold without boundary, and 7 a simple 



non-degenerate geodesic on M, with the spectrum of its linear Poincare map satisfying (2.9). 
Then the spectrum of the Laplacian, Ag, determines the Birkhoff normal form of"f. 

In the special case of symmetric domains in the plane, we can apply Theorem ^ and an 
observation of Colin de Verdiere (H §4] to recover a result of Zelditch pi|p|: 

Corollary 5.5. Suppose that O C is a bounded domain, dfl is real analytic, and that fl is 
symmetric with respect to the x and y axes. Suppose also that has a simple non- degenerate 
bouncing ball orbit along one of the axes. Then the spectrum of the Dirichlet ( or Neumann ) 
Laplacian determines f2 among all such domains. 



^Alternatively, we could have reduced the problem to the boundary, and considered the monodromy operator 
there. 

^Melrose |^] has announced the same result, sketching a yet different proof. Recently, by considering the full 
quantum Birkhoff normal form, Zelditch succeeded in removing one of the symmetry conditions 
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